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Abstract 

The dynamics of the Reynolds stress tensor is determined by an evo- 
lution equation coupling geometrical effects and turbulent source terms. 
The effects of the mean ffow geometry are shown up when the source 
terms are neglected. Then, the Reynolds stress tensor is expressed as the 
sum of three tensor products of vector fields, which are governed by a dis- 
toTted gyroscopic equation. Along the mean flow trajectories and in the 
directions of the vector flelds, the fluctuations of velocity are determined 
by differential equations whose coefficients depend only on the mean flow 
deformation. 



1 Introduction 

The governing equations of barotropic turbulent compressible fluids are (see, for 
example [J, g], [3]) 

{P)t + {{P) U^)n=0, 
Hp) U^)t + ((P) U^UJ + (p) % + {pu,Uj)) „ = 0, (1) 

{puiU.j)^ + {{puiUj) Uk) ,k + (pUiUj) U-j,k + (pu^Uj) Ui^k = S^j . 

Here "brackets" mean the averaging, "coma" means the derivation with respect 
to the Eulerian coordinates x = {x^}, i G {1,2,3}, index t means the partial 
derivative with respect to time, p is the fluid density, U = {f/i}, i 6 {1, 2, 3} is 
the mass average velocity, p is the pressure, u = {ui}, i G {1, 2, 3} is the velocity 
fluctuation verifying (pu) = 0. Repeated indices mean summation. The term 
S = {Sij} represents turbulent sources. We introduce the Reynolds stress tensor 

R=(/9U(g)u), [Rij = {puiUj)). 
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The system ((TJ can also be rewritten in the tensorial form 



d{p) 
dt 



div((p) U) = 0, 



(p)^ + V(p) + (divRf 



dR „ , _ dlJ ^ ^ fdV 
-+Rd.vU+-R + R - 



(2) 



where d/dt means the material derivative with respect to the mean motion 

d d rp 

The superscript " " means the transposition. By using the mass conservation 
law, the equation for the volumic Reynolds tensor R can be rewritten as the 
equation for the specific (or per unit mass) Reynolds tensor 



dt dx \dxj (p) ' 



(3) 



where 



The aim of this note is to understand the structure of homogeneous Reynolds 
equation ([3]) corresponding to S = 0. The reason is twofold: 

- first, in the numerical study of compressible turbulent flows, this is a natural 
step in applying the splitting- up technique (see for example [1]), 

- second, the system ([5]) also appears as an exact asymptotic model of weakly 
shearing flows of long waves over a fiat bottom [5] 



dh 
'dt 



+ div( hU) = 0, 



dt V 2 



(divR)^ = 0, 



(4) 



— +RdivU+ — R + R — 
at ox \ ox 



0. 



In Eq. h is the fluid depth, the average pressure is given by (p) — gh? g 
is the gravity acceleration, and 



R 



(U - U) (8) (U - \])dz, h\] ■ 



15 dz, 
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where U is an instantaneous velocity. Equations are written for three-dimensional 
long waves and the production term is zero in the limit of weakly shearing flows. 
Also, Eqs. (H]) are hyperbolic. Finally, we will focus on the equation of the 
Reynolds tensor per unit mass 



The particular case rot U = was investigated in [6]. In such a case (9U/9x)^ = 
d\J/dx and Eq. ([5|) corresponds to a two-covariant tensor convected by the 
mean flow. This means that P has a zero Lie derivative with respect to the 
velocity field U and the tensor Pq, image of P in Lagrange coordinates {t, X), 
only depends on X = {Xi}, i G {1, 2, 3} 

rf.P ^ f + 1^ P + = 0, P = {F-y' Po (X) F'\ 

dt ox ox ^ ' 

where F = i9x/9X is the deformation gradient of the mean motion [7]. 

The aim of the paper is the study of the homogeneous Reynolds tensor 
equations structure ([5]) in the case rot JJ ^ 0. 



2 Geometrical properties of the Reynolds tensor 
evolution 

The Reynolds stress tensor P is symmetric, semi-positive definite and conse- 
quently can be rewritten in a local basis of orthonormal eigenvectors in the 
form 

3 3 
a—1 a—1 

The eigenvalues A^, a € {1, 2, 3} are non-negative; the case A^ > is a generic 
one. For the two-dimensional case, A3 = 0. Let us denote 

aa = Xa Ba, (Aq > 0) « G {1, 2, 3}. 

Then, 

3 3 
P ^ aa (g) Ha = ^ aa a^. (6) 
a—1 a=l 

From Eq. jB]), we deduce 

a—1 ^ 

By using Eq. ((7|), Eq. ^ can be written 



The vector daa / dt + {dXJ / dx.) can be developpcd in the local basis {aa}, en G 
{1, 2, 3} of eigenvectors; one obtains 

dan dXJ . , , 

^ + ^ a„ = }J^;3ca^, a e {1,2,3} (9) 

^ /3=1 

where A/^a, (a,/? 6 {1,2,3}) are the scalar components to be determined. By 
using Eq. ©, Eq. ® leads to 

3 3 
^^QQgQQ+ ^ {Aap + Apa) gQ/3 = 0, 
Q=l a#/3=l 

where g„c< = 23^ and g^^a = g^j^ = a^ aj + a^, a^, (a, /3 G {1, 2, 3}), are 
six independent symmetric tensors. Consequently, 

Aaa^O and ^c«/3 + ^/3q = 0, a, /3 e {1,2,3}. 

Equation ([5]) is equivalent to 

^ + ^ a„ = Ai (tt) with tt =^32 ei + A13 e2 + ^21 63, a € {1,2,3}, 
dt ox 

. 

where diagonal matrix A and antisymmetric matrix i (tt) are determined in the 
basis Gq, a £ {1, 2, 3} as 

Ai 
A = I A2 
A3 

The vectors a^,, a G {1, 2, 3} are orthogonal, a^ a^ =0, (a / If we assume 
the orthogonality at initial instant, this is equivalent to 







-^21 


^13 


^(7^)= 1 


^21 





-A32 




V -^13 


^32 






So, Eqs. (Uni) - dni) yleld 



Va7^/?G{1,2,3}, 



a„ ( Ai (7r)e/5 - — a^ 1 +a^ I A z (tt) - a^ 



or 

2Aa A/3eJ D a^ A i (tt) e/3 + a^ A i (tt) 

(12) 

= A^e^ i (tt) + A^e^ i (tt) e„ 



4 



where 

D 



1 / 9U /OU^ ^'^ 



2 \ dx \d: 



is the rate of deformation tensor corresponding to the mean flow. We denote 
the mixed product of three vectors {a, b,c} as (a, b,c) = a^(bAc). Hence, 
Eq. ini) is 

2XaXf3 e^B ep = A^) (e„, tt, e^) 

= ( -^^-^a) (7i-,e„,e^) 

where {a, 7} is a cycUc permutation of the triplet {1, 2, 3}. Finally, we get 

-^e, . ^^"^^r\?^^ (13) 
Equation (|10p can be written 

+ + = ae{l,2,3}. (14) 

at at ox 



Since 



^ = 0, e^A i (tt) = A^e^i (tt) = 0, 



by multiplying the left side of Eq. with e^, we get 

^+ (elDe„)A„ = 0. 

By multiplying the left side of Eq. p4|) with the projector (I — eae^), we get 

dea , d\J , . / X 

Aq— + Aq (I - Gae^) — -Ai(7r)ea=0. 

Due to the fact there exists a vector 11 such that 

den 



vector n verifies the condition 



— = n Ae„, 
dt 



Xai (U) Ba + Xa {I ~ eael) — Ba -Ai(7r)eQ = 0. (15) 
By multiplying Eq. (jlSp with where /3 7^ a, we get 

Aall'^e^ = A^Tr'^e-y - -^"^^"^ (l^) 
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By replacing Rel. (fT^ into Eq. (fTC)) we get 



2„T 



ax v 9x 



e/3 



— 



We can now formulate the following result: 

Theorem 1. The Reynolds tensor can be written in the form 

3 



R = (p) ^ ea ® ea. 



The eigenvectors and the eigenvalues Aq, verify the equations: 
den 



< 



dt 



= n A gq, 



d (Ln Aq 
df 



-2 Ma> 



a £{1,2,3} (17) 



where 



Ma = D e„, n^e^ 



+ Xi— I ea 



A/3 - A„ 



The triplet {a, (3, 7} corresponds to a cyclic permutation of the triplet {1, 2, 3} . 

Equation pTl *) is similar to the equations of a rigid body [8J. The vectors 
Gq, form a natural moving frame G = {^0)^1=1 depending on the mean deforma- 
tion. The eigenvalues A^ of the Reynolds tensor are determined by the evolution 
equation (|17b ). Note that if Aq are initially positive, they will be positive at 
any time. Hence, the tensor P is always positive definite. 

Due to the mass conservation law ([21) and Eq. (fTTl ). we obtain 
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Consequently, system ([T7|) admits a scalar invariant along the trajectories of 
mean flow. This invariant was earlier obtained in f6] in a different form. 

Let us introduce the turbulent specific energy 

3 



In the incompressible (isochoric) case, we have d{p) jdt — 0; the turbulent en- 
ergy is minimal in the isotropic case when the three eigenvalues A„ are equal 
(Aj = A2 = A3 = A ). 

In this case, the orthonormal eigenvectors e^, a £ {1,2,3} of the Reynolds 
stress tensor P are also the orthonormal eigenvectors of the mean rate of defor- 
mation tensor D. The eigenvalues of D are /i^, a S {1,2,3}. 

2 /3 2 2 /3 

In the compressible isotropic case ct = (p) k, k = 3A /(2(p) ' ), and 
K is a classical invariant of isotropic turbulence. In presence of shock waves 

-2 ^ 

the quantity (p) Yl iio^ conserved through shocks; it increases like the 

a=l 

classical entropy in compressible fluid dynamics. The estimation of the jump of 

turbulence entropy in isotropic case was given in [9]. 

The governing equations ([2]) admit the energy conservation law 

I ((p) Q|U|^ + e, + er) ) + div (^(p)U Q |U|^ + e,) + + R) = 0, 

where the internal specific energy is defined by 

1 



de^ = -{p)d ^ 

and the mean pressure (p) is supposed to be a given function of (p). Indeed, 
using pTl ) we immediately obtain 

I ((p) Q|U|^ + + er)) + div (^(p)U Q|U|^ + e,) + {{p}! + R) U 

(p)^+tr(RD) = 
(p)tr 





3 Conclusion 

The equations of fluid turbulent motions take three equations into account: the 
equation of the mass balance ([H), the balance equation of the average momen- 
tum (IH), and the Reynolds stress tensor equation of evolution ([T|); this last 
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equation has been the object of our study. If the turbulent sources are neglected, 
the turbulent fluid motion is a superposition of the mean motion and turbulent 
fluctuations. The eigenvectors of the Reynolds tensor carry the fluctuations 
associated with the mean flow deformation. The amplitude of turbulent defor- 
mations is defined by the eigenvalues of the Reynolds stress tensor. Locally, 
the equations for the directions of turbulent fluctuations describe a small solid 
whose rotation is given by a gyroscopic type equation (Equation pTl )). The 
amplitude evolution of turbulent deformations is determined by the diagonal 
values of the mean deformation tensor D expressed in the eigenvector basis 
of the Reynolds stress tensor. The turbulence increases with the time when 
< 0, and decreases when > 0. In the particular case of incompress- 
ible fluid motions we have tr D = 0, and hence there always exists a direction 
in which the turbulence is increasing while in other directions it is decreasing. 
These mathematical deductions are confirmed by experiments |lOj . 
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